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Abstract. We study the initial-boundary value problem of the Navier-Stokes equations 
for incompressible fluids in a general domain in R n with compact and smooth boundary, 
subject to the kinematic and vorticity boundary conditions on the non-flat boundary. 
We observe that, under the nonhomogeneous boundary conditions, the pressure p can 
be still recovered by solving the Neumann problem for the Poisson equation. Then we 
establish the well-posedness of the unsteady Stokes equations and employ the solution 
to reduce our initial-boundary value problem into an initial-boundary value problem 
with absolute boundary conditions. Based on this, we first establish the well-posedness 
for an appropriate local linearized problem with the absolute boundary conditions and 
the initial condition (without the incompressibility condition), which establishes a ve- 
locity mapping. Then we develop apriori estimates for the velocity mapping, especially 
involving the Sobolev norm for the time-derivative of the mapping to deal with the 
complicated boundary conditions, which leads to the existence of the fixed point of the 
mapping and the existence of solutions to our initial-boundary value problem. Finally, 
we establish that, when the viscosity coefficient tends zero, the strong solutions of the 
initial-boundary value problem in R n (n > 3) with nonhomogeneous vorticity boundary 
condition converges in L 2 to the corresponding Euler equations satisfying the kinematic 
condition. 



1. Introduction 

The motion of an incompressible viscous fluid in M n , n > 2, is described by the Navier- 
Stokes equations: 

dtu + u ■ Vu + Vp = fiAu, (1.1) 

V-ti = 0, (1.2) 

with initial data 

u\ t =o = u (x), (1.3) 

where u(t,x) = («,••• ,u n )(t,x) is the velocity vector field and p(t,x) is the pressure 

that maintains the incompressibility of a fluid at (t,x). Equation (jl.2p . i.e. div u = 0, 
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is the incompressibility condition. As a nonlinear system of partial differential equations, 
u and p are regarded as unknown functions, and the initial velocity field uq(x) sets the 
fluid in motion. The constant n > is the kinematic viscosity constant, u ■ V denotes the 
covariant derivative along the flow trajectories, namely, the directional derivative in the 
direction u, Au is the usual Laplacian on u, and fiAu represents the stress applied to the 
fluid. As usual, we use V- = div to denote the divergence operator. 

For inviscid flow, u = and then the equations are referred to as the Euler equations 
for incompressible fluid flow: 

d t u + u ■ Vu + Vp = 0, (1.4) 
V-u = 0. (1.5) 

When a fluid is confined in a bounded domain £1 C M. n with non-empty boundary T, 
these equations must be supplied with proper boundary conditions in order to be well- 
posed. For concreteness, the bounded domain Q C R n is assumed to have a compact, 
oriented, smooth surface boundary T. In this paper, we propose and study the initial- 
boundary value problem (jl.lj) - (jl.3j) with the following boundary conditions on the non-flat 
boundary F: 



u 



V 



(kinematic condition), (1-6) 
a (vorticity condition), (1-7) 



where the field a = a(t, x) is defined on the non-flat boundary r, S R denotes the 
normal component and the tangential part of u £ W 1 . 

Various physical considerations for some flows observed in the nature have led to the 
no-slip condition: 

nil = 0, (1.8) 

together with the kinematic condition fjl .6[) . This is equivalent to the simple homoge- 
neous Dirichlet boundary condition: u|r = 0. Thus, it has been received an intensive 
investigation for the Navier-Stokes equations. 

On the other hand, if a vector field u £ M. n is identified with its corresponding dif- 
ferential 1-form on Q, then various boundary problems for the linearized equations, the 
heat equations, have been studied in geometry under the Hodge theory on manifolds with 
boundaries; see Conner [9] and the references cited therein. In particular, three natural 
boundary conditions have been identified; since 

Au = -Vxw + V(divit) onQcK", 

where u> := V x u £ M. n is the vorticity, boundary conditions may be imposed on the 
components u^, n" , u^, w" , and div u. It has been revealed that the initial-boundary value 
problem with the boundary conditions (jl.6p ~ (jl.7p for the heat equations is well-posed. 
This leads to a natural question whether the same setup of initial-boundary value problem 
is well-posed for the Navier-Stokes equations. In the literature, the absolute boundary 
conditions mean that our boundary conditions (11.6l )-( fl~7|) with a = 0, i.e. homogeneous. 

The initial-boundary value problem (|l.lj) - ([1.3|l and (|1.6j) - (jl.7|l for the Navier-Stokes 
equations has been received less attention, partly due to the common agreement that the 
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no-slip condition is a physical condition under ideal situations. However, recent experi- 
mental evidence (cf. Einzel-Panzer-Liu [11], Lauga-Brenner-Stone [23], and Zhu-Garnick 
|55j for a survey) demonstrates that the velocity field of a fluid does not satisfy the com- 
plete no-slip conditions in general. Indeed, it is revealed that the pressure does depend on 
how curved the surface boundary is; also see Bellout-Neustpua-Penel [6] . It is also evident 
that, for high speed flows, one can not expect to comply the complete no-slip conditions. 
Therefore, a careful study of incompressible fluids with various boundary conditions asso- 
ciated with the Navier-Stokes equations is in great demand. 

Another motivation for our study in this paper is from the Navier slip boundary con- 
dition, which are widely accepted in many applications and numerical studies. It states 
that the slip velocity is proportional to the shear stress: 

u ll =-(((Vu + (Vu) T )v)\ (1.9) 

where v is the outward normal on Y and C > is the slip length on T. Such boundary 
conditions can be induced by effects of free capillary boundaries, a rough boundary, a 
perforated boundary, or an exterior electric field (cf. [H IU EJ [TJ EJ H5J [HJ S3]). Take an 
asymptotic expansion of u(t,x) with respect to ( and the distance d(x,T) to T such that 
d(x, r) <C Then, by a direct calculation, one finds that the Navier slip condition (|1.9j) 
yields (|1.7p in which a(t, x) is determined by the initial value problem of a family of linear 
heat equations on T, that is, a(t,x) is determined solely by the initial data fjl .31) . 

There is another appealing reason to consider other boundary problems associated with 
the Navier-Stokes equations. The Euler equations (|1.4|) - (|1.5p subject to the kinematic 
condition (|1.6p , along with the initial condition (|1.3p , are often regarded as an ideal model 
in turbulence theory, where the Reynolds number is so big that the viscosity coefficient 
fi is negligible. Hence, a natural question is whether the flow determined by the Euler 
equations is, at best, a "singular limit" of the Navier-Stokes equations as the viscosity 
coefficient \i \ (cf. Pope [39], page 18). However, it is very subtle and difficult in general, 
since physical boundary layers may be present. The classical no-slip boundary condition, 
u\y = 0, give arise to the phenomenon of strong boundary layers in general as formally 
derived by Prandtl [40] . For some sufficient conditions to ensure the convergence of viscous 
solutions to the ones of the Euler equations, see Kato [IHl [52] and the references therein. 
However, as far as we know, such a claim has never been proved for general bounded 
domains with the Navier boundary condition. In contrary, we prove in this paper that 
our problem (|1.1I) - (|1.3I) and (jl.6p - (jl.7p for the Navier-Stokes equations converges in L 2 
to a solution of problem (|1.4p ~ (|1.6p for the Euler equations as the viscosity coefficient 
(j, I 0. To explain why the Euler equations, subject to the kinematic condition (11. 6h . is 
the limiting case of the solutions u M of our problem (|l.ip - (jl.3p and (|1.6p - (|1.7p . we notice 
that the pressure p required in the Euler equations solves the Poisson equation: 

Ap = — V • (u ■ Vu) , d u p\ r = it(u,u), (1-10) 

where u(t, x) is a solution to the Euler equations (|1.4p - (|1.5p . and tt is the second funda- 
mental form of the boundary T which describes the curvature of the boundary surface 
r. Since u does not necessarily comply the no-slip condition, the normal derivative of 
the pressure p does not vanish along the boundary in general. In particular, the pres- 
sure p along the boundary naturally depends on the curvature of the boundary. On the 
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other hand, under the kinematic condition (jl.6p the pressure p can be recovered from the 
Neumann problem of the Poisson equation: 

Ap = -V-(« M -V^), (1.11) 
d»p\r = niu^Ufj) - fiV r x wjl. (1.12) 

As n [ 0, the pressure p satisfies the same equation as that of the pressure required by 
the Euler equations. However, if we further enforce the no-slip condition, then the normal 
derivative of the pressure p vanishes, which does not in general coincide with that of the 
Euler equations. On the other hand, if we set a free condition on the tangent component 
of u but replace the boundary condition on the vorticity, we have the chance to recover 
the right boundary condition for the pressure in the Euler equations. 

The rigorous mathematical analysis of the Navier-Stokes equations involving vorticity- 
type boundary conditions may date back the work by Solonnikov-Scadilov [36] for the 
stationary linearized Navier-Stokes system under boundary conditions: 



u 



0, ((Vu + (Vu) T )z/)"| r = 0. (1.13) 



For two-dimensional, simply connected, bounded domains, the vanishing viscosity problem 
has been rigorously justified by Yodovich [33]. We also refer to Lions [281 [29] for the 
boundary conditions (jl.6p - (jl.7p with a = 0, Clopeau, Mikelic, and Robert [8] and Lopes 
Filho, Nussenzveig Lopes and Planas [30] for the Navier boundary condition (|1.13p . and 
Mucha [35] under some geometrical constraints on the shape of the domains for the two- 
dimensional case. Also see |53j for an analysis for the complete slip boundary conditions 
as (PJ3D below. 

The main purpose of this paper is to develop an approach to deal with the well-posedness 
and the inviscid limit for the initial-boundary value problem (|l.lj) - ([1.3j) and (jl.6f> (|1.7h 
with nonhomogeneous boundary condition for the multidimensional Navier-Stokes equa- 
tions (n > 3). One of the main difficulties for the higher dimension case, in comparison 
with the two-dimensional case, is that the maximum principle for the vorticity fails so that 
the two-dimensional techniques can not be directly extended to this case. Furthermore, the 
nonhomogeneous boundary condition causes another main difficulty in developing apriori 
estimates which require to be compatible with the nonlinear convection term. In Section 
2, we introduce the local moving frames on the boundary T of a domain O and derive 
some basic identities for vector fields. Although we work with a domain in the Euclidean 
space, the boundary T is a curved surface, so that geometric tools have been brought in 
to carry out local computations. Since we can work only with the boundary coordinate 
system which may be not "normal" , the computations we need to carry out are a little bit 
complicated and long. The results in Section 2 allow us to determine boundary values of 
several interesting quantities, which will be used to settle the local existence of a unique 
strong solution for problem (jl.ip - (jl.3p and (|1.6|) - ([1.7|) for the Navier-Stokes equations. 

In Section 3, we establish various L 2 -estimate for vector fields that satisfy the absolute 
boundary conditions: 

u 1 - = 0, J = 0, (1.14) 
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In particular, the main novel fact is that, if u(t, x) satisfies the absolute boundary condi- 
tions (|1.14p . then 



along the boundary is crucial in obtaining the necessary apriori estimates for the lineariza- 
tion of the absolute boundary problem. 

In Sections 4, we develop an old idea of Leray-Hopf [HI [24J, [25j [26] for our initial- 
boundary value problem with nonhomogeneous vorticity boundary condition for the Navier- 
Stokes equations. Namely, under the nonhomogeneous boundary conditions, the pressure 
p may be recovered by solving the Neumann problem (jl.lip - (|1.12p of the Poisson equa- 
tions. Then we establish the well-posedness of the unsteady Stokes equations and employ 
the solution w to reduce our initial-boundary value problem into an initial-boundary value 
problem with absolute boundary conditions for v = u — w. This leads to the following 
linearlization: One first considers the linear parabolic equation: 



where pp solves ([l.lOp replacing u by (3 + w for a given f3 = f3(t, x),t < T, satisfying the 
absolute conditions and the initial condition. Based on this, in Section 5, we demonstrate 
that the absolute boundary conditions are suitable boundary conditions, since they gen- 
erally do not appear in the literature on parabolic equations, and conclude that (|1.15p 
is well-posed, which thus establishes a velocity mapping V : w — ► v. Note that we have 
dropped the incompressibility condition, so that (|1.15p is local and a linear parabolic 
equation. However, we will establish that any fixed point in a proper functional space 
is a strong solution to the Navier-Stokes equations, so that the incompressibility can be 
recovered. 

We emphasize that the absolute boundary conditions in a general domain are a kind of 
boundary conditions which change from the Dirichlet to the Neumann boundary condition 
along the boundary. Therefore, the coercive estimates, or called the global L p -estimates, 
for parabolic equations are not applicable in this setting. Indeed, we have to establish 
estimates quite close to the coercive estimates for the linear systems (j!.15j) in a different 
functional space, which will be done in Section 6. In order to construct a local (in time) 
strong solution, we need to develop apriori estimates for the solution v(t,x) of the linear 
parabolic equation (j!.15|) . It would be enough to develop an estimate for the i? 2 -norm 
of v(t,x) in terms of that of (/? + w)(t,x). In fact, it works for the Dirichlet boundary 
problem (which has been done by Solonnikov [45] ) and also works for the periodic case 
(cf. [19J). Unfortunately, due to the complicated boundary conditions, we are unable 
to achieve such apriori estimates, instead, we have to involve the Sobolev norm for the 
time-derivative of v(t,x). 

Finally, in Sections 8-9, we establish that the unique strong solution of the initial- 
boundary value problem f|l.l[) — f)1.3j> and (|1.6I) - ([1.7I) in M n ,n > 3, with nonhomogeneous 
vorticity boundary condition converges in L 2 to the corresponding Euler equations satis- 
fying the kinematic condition (II. 6p . 



(V x u) 1 - = 0. 



The latter information on 



V x oj = Au — V(div u) 




(1.15) 
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2. Local Moving Frames on the Boundaries and Basic Identities for 

Vector Fields 

Let Q C R n be a bounded domain with a compact, smooth boundary T = d£l with 
n > 2. We assume that T is an oriented surface (not necessarily connected) which carries 
an induced metric. Let v denote the unit normal vector field along T pointing outwards. As 
usual, we use conventional notation that the repeated indices in a formula are understood 
to be summed up from 1 to n unless confusion may occur. 

In order to reflect better the geometry of the boundary T = d£l of f2, it is convenient to 
work with local moving frames compatible to T (see Palais- Terng [37j for the details) . More 
precisely, by a moving frame compatible to the boundary we mean a local orthonormal 
basis 

{ei,--- ,e n } 

of the tangent bundle TQ of Q about a boundary point x 6 T such that e n = v when 
restricted to T, where v is the unit normal vector field along T pointing outwards. 

The Christoffel symbols are defined by Y l - = OUei). In terms of the Christoffel symbols: 

V et ej = r^ez, V ez a j = -Y^a 1 , 

where Vj means the directional derivative in X. The torsion- free condition may be stated 
as T^. = — T J ik and, along the boundary T = dQ, 

T™- = = — hij, for any 1 < i,j < n — 1, 

and 7r = {hij) is the second fundamental form, which is a symmetric tensor on T. By 
definition, 

hij = {V e M,ej). 

If u = ^27=1 u%e i an d w = YllZi w%& i are two vector fields tangent to T, then 

n-1 

it(u,w) = hijtfwi . 

The surface T carries an induced metric, its Levi-Civita connection is denoted by V r . 
Then {ei,-- - ,e n _i} restricted to T is a moving frame of TT. If u is a vector field on 
ft, then the tangent part of u restricted to T, denoted by it II, is a section of TT, and its 
covariant derivative is denoted by V r n". Let H denote the mean curvature of T, that is, 
H is the trace of the second fundamental form n: 

n-1 
3=1 

Lemma 2.1. Let u and w be two vector fields on Q. Then 

— 7r(u",u;") — H(w, v){u, v) + (w, v)(V • u) 
+(J,V r (it,i/)) + (J ,V r (w,u)) -V r - ((w, u)J). (2.1) 
= on V, then 

(w-Vu,v)=—k(u,w). (2.2) 



(w ■ Vu, v) = 
In particular, if u 1 - = w 1 - 
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Proof. Since 

{wVu) l = wi{e 3 {v?)+T) k u k ), 

then 

(w ■ Vu, v) = w j ( ej (u n ) + T] k u k ) 

= (w\V T (u, v)) -ir(J, W W) + w n {V n u n ). (2.3) 

On the other hand, according to the Ricci equation, 
n— 1 Ti—1 Ti—1 

\7 jU j = Vj« j + ^ /^(u, i/) = V r • ull + H(u, u), 
i=i i=i i=i 

so that 

n— 1 

V n u n = V • it - ^ Vyu* = V • u - V r • u 11 - H(u, v). (2.4) 
i=i 

Substitution (|23J) into ([13]) yields 

(w • Vit, i/) = (to 11 , V r (u, i/)) - ^(u 11 ,^ 11 ) + (w,u)(V • it - V r • u 11 - iJ(w, v)), 
which, together with the identity: 

{w,v)V r ■ J = V r • ({w,v)J) - {J,V r {w,v)), 
yields ([23]) . □ 

Lemma 2.2. Lei u be a vector field on Q, lo = V x u, and d = V • u. Then 

(i) c^cilp = (Ait, zv) + V r x cjll, u;/iere 

V r x J 1 = V\uj 2 - V\uj 1 

is independent of the choice of a local moving frame, which can be identified with the 
exterior derivative o/w" on T; 

(ii) |di/(M 2 )|r = {u x. lo, v) + {u, v)d — 7r(?J, «H) — # (it, z^) 2 + 2(JI, V r (-u, v)) 

-V r • ((u, i/)itH). 

Proof, (i) According to the vector identity: 

Au = -V x lo + Vd, 

one obtains 

cWlr = (Au, i/) + (V x lo, v). 

It remains to verify that (V x lo, v) coincides with V r x which depends only on the 
tangent part cJ'| r . The last fact follows easily from a local computation, 
(ii) This formula follows directly from the vector identity: 

I\7(|n| 2 ) = u x (V x u) + (u ■ V) u, 
and Lemma l2.1i □ 
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3. l 2 -estimates for vector fields satisfying the absolute boundary 

Conditions 

In this section, we establish L 2 -estimates for vector fields satisfying the absolute bound- 
ary conditions (|1.14p . 

Let C I 3 be a bounded domain with a compact, smooth boundary T. The Sobolev 
spaces can be defined in terms of the total derivative V. If T is a (smooth) vector field on 
O and V J T is the j-th covariant derivative, then, for a nonnegative integer k, the Sobolev 
norm || • || W fe, P ,p > 1, is given by 

\\n w ^ = {jz\\^n%) 1/p , (3.i) 

3=0 

where || • \\lp denotes the L p -norm of vector fields, that is, = (J*^ \T\ p dx) 1 l' p with 

\T\ = y/{T, T) the length of the vector field. When p = 2, we often use H k (n) = W k > 2 {Q) 
and ||T||^fc = ||T|| W fc,2, since such spaces are Hilbert spaces. The most important in the 
theory of Sobolev spaces is the Sobolev imbedding theorem (see Theorem 4.12 of Adams- 
Fournier j2J 3 page 85). 

We will use the following integration by parts formula: 



J (u x w,u)dS, (3.2) 



(V x u, w) dx = I (u, V x w) dx + / (u X w,v) dS. 
Jn Jr 

where dS is the surface measure on T. This formula follows easily from the Stokes theorem. 
If u is a vector field on Q, then three L 2 -norms are related to the total derivative V«: 

(i) The (squared) L 2 -norm HViilH; 

(ii) ||V-u|| 2 + ||Vxu||2, which is most useful in the study of the Navier-Stokes equations, 
since, if u is the velocity field of a incompressible fluid, then V • u = and V x u = to is 
the vorticity in physics; 

(iii) The quantity — J Q (Au, u) which is related to the spectral gap of the (Hodge De 
Rham) Laplacian A which has a clear geometric meaning. 

The relations among these three quantities are the following. 

Lemma 3.1. If u £ H 2 (£l) is a vector field, then 

/ {Au,u)dx = — I \uj\ 2 dx — / |V • u\ 2 dx 
Jn Jn Jn 



and 



+ j (u x uj,v)dS + j (V -u)(u,u)dS, (3.3) 

/ |Vn| 2 dx = / |u| 2 dx+ / |V • u\ 2 dx — / 7r(u", u") dS 
Jn Jn Jn Jr 

-J H\u ± \ 2 dS + 2 J (J,V r (u,u))dS. (3.4) 

Proof. Formula (|3.3|) follows from the vector identity 

Am = -Vxw + V(V-u) (3.5) 
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and integration by parts formula (|3.2p . To show (|3.4p . integrating the Bochner's identity: 

(Ait.it) = ^A|n| 2 - |Vu| 2 , (3.6) 
together with integration by parts formulas, we obtain 

f (Au,u) = - [ \Vu\ 2 dx+- [ d v {\u\ 2 )dS 
Jn Jn 2 Jy 

= - |Vu| 2 dx + / (n 11 x <J,v)dS + (u,v)V-udS 
Jn Jr Jy 

+2 J (J,V r (u,v))dS - J ir(J,J)dS-J H(u,v) 2 dS, (3.7) 



where we have used Lemma l2.2f ii). Now (|3.4p follows from (|3.3p and (|3.7p . □ 

Lemma 3.2. Suppose that u G H 1 ^) is a vector field satisfying the boundary conditions 
that u"L = and uM„ = 0, where T = T± U To. Then 

III II 2 7 

H ff i<C||(w,V-«,«)|| a (3.8) 
/or some constant C depending only on fi. 

Proof. According to (|3.4p and the boundary conditions, we have 



/ |Vu| 2 dx = / \(u,V ■ u)\ 2 dx - [ ir(J,J)dS- [ H(u,v) 2 dS 
Jn Jn Jr 2 Jtx 



< / -uTdx + C / \u\ 2 dS. 



Now the conclusion follows from the trace imbedding theorem (Theorem 1.5.1.10 in Gris- 
vard [12J): For any e £ (0, 1], 

\u\ 2 dS<e [ \Vu\ 2 dx + — I \u\ 2 dx, (3.9) 



e Jn 

where C is a constant depending only on 0. □ 

Remark 3.1. The above lemma holds not only for the L 2 -norm, but also for the L p -norm 
for any p S [2,oo). Indeed, under the same boundary conditions on u as in Lemma 13.21 
as a special case of Theorem 10.5 of [3], we have 

\\u\\ w i, P < C(||(w,V-u)||p + H| 2 ), (3.10) 

where C also depends on p > 2. Estimate ()3.10p is the Agmon-Douglis-Nirenberg's esti- 
mate (cf. [3]). 

Next we consider the second derivative V 2 u of u. 

Lemma 3.3. There exists C > depending only on £1 such that 

\\V 2 u\\l < C(||Au||| + ||V«|||) 

for any vector field u £ H 2 satisfying the absolute boundary conditions (|1.14p . 
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Proof. Apply Lemma l2.2l f ii) to a gradient field V/ for any scalar function / on £1. Then 

~d,(|V/| 2 ) = -tt((V/)II, (V/)H) - H \d w f\ 2 + d„fAf 

+2((V/)H , V r (d u f)) - V r • {d v f (V/)" ). (3.11) 

Write u = (u , • • • , u n ) under an orthonormal basis of TQ. Then, according to the Bochner 
identity (pTBjh 



n n i 

|V 2 n| 2 = |Wp = ^ (^A| Wf - (VA/,V/)), 

fc=l k=l 

and, after integration, we find 

\\V 2 u\\l = Vf- / A|Vn fc | 2 dx- / (VA^.V/)) efe 

= J2 I (Au k ) 2 dx+- [ d u (\Vu k \ 2 )dS - [ (Au k )d u u k dS 
j^Jn 2 Jr Jr 

n ,. ,. n—1 

= V / {Au k ) 2 dx- / S2 hijiViU^iVju^dS 
k=i Jn Jr i,j=l 

- [ Y\d u u k \ 2 HdS + 2 [ (V r u k ,V r (d u u k ))dS 

/n— 1 „ n 

h ij V i u k V j u k dS - / ^2(Vu k ,u) 2 HdS 

' i,j=l r k=l 

n „ 

+2V / (V r u k ,V r (Vu k ,u))dS + \\Au\\l 

where the third equality follows from (|3.1ip applying to each u k . We now handle the 
boundary integral: 

V [ {v r u k ,v r (d u u k ))ds = V / {v r u k ,v r (d u u k ))ds + I {v r u n ,v r (d u u n ))ds 
k=l Jr k=l Jr Jr 

where we have used the fact that u"|r = so that S7 T u n = 0. For 1 < k < n — 1, we have 
V n u fc = VfcM™ = on T. However, 

n—1 n—1 

nj> 



du u k = V n U k - ^nj = ~ E 



3=1 3=1 
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so that 

n— 1 n— 1 n— 1 

£<v r «* v r (<*,u*)) = -£< vr « fc . v r ^>i%.) < c (H 2 + |Vn| 2 ) , 
fc=i fc=i j=i 

and the inequality follows the trace imbedding theorem. □ 

Corollary 3.1. Let u G H 2 (Q) be a vector field which satisfies the absolute boundary 
condition (|1.14p on T. Let d = V • u, w = V x u, and i\) = V x a;. TTien 

(i) V • w = V • V = 0; 

(ii) An = Vd - V x 

(hi) u ± \ r = 0, w"| r = 0, and V^lr = 0; 

(iv) There exists C > depending only on f2 suc/t t/iai 

||«||j?a < C||(Vd,V X uj,u)\\ 2 . 

That is, ||(Vd, V x cj,u)||2 is an equivalent norm for any vector field u G id 2 (f2) satisfying 
the absolute boundary condition (jl,14p . 

Identity (i) is obvious, and (ii) follows from the vector identity: 

Au = - V x (V x u) + V (V • u) . 
To show V^lr = 0) we use (ii) to obtain 

(Au) ± = (Vd) ± -i; ± , 

and the claim follows from Lemma 12,21 According to Lemma 13.31 and the Ehrling- 
Nirenberg-Gagliardo interpolation inequality, one has 

\\u\\ H 2 < C(||Au|| 2 + \\uW2), 

and (iv) follows from (ii). 

4. Navier-Stokes Equations 

Consider the Navier-Stokes equations (jl.ip - (jl.2p in a bounded domain f2 C M 3 , subject 
to the kinematic condition (jl.6p and the nonhomogeneous vorticity condition (jl.7p . 
Using the vector identity: 

u- Vu = uj x u + ^V(|u| 2 ), (4.1) 

one may rewrite (jl.ip as 

d t u + u x u = fiAu - V(p + ^H 2 )- (4.2) 
4.1. Neumann problem for the pressure p. 

Lemma 4.1. Ifu(t,x) is a smooth solution of the Navier-Stokes equations (|l.ll )-( Pl~2|) m 
$7^ := x (0, T] subject to the boundary conditions (|1.6j) - (|1.7p . thenp solves the Neumann 
problem: 

Ap = - V • (u • Vu) , d„p| r = vr(u, u) - uV r x a. (4.3) 



12 
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Proof. Since u(t,x) is a smooth solution in Qt, by taking the normal part of equations 
(jl.ip . we may determine the normal derivative of the pressure p. Indeed, since u ± \ r = 0, 
equations (jl.ip imply that 

((u • V)u) ± = (Au) — (Vp) ± on the boundary T. 

We apply Lemmas I2.1H2.2I to obtain 

9uP\r = p{Au, v) — (u- Vu, u) 

= pd u (V ■ u) — p(V x uj, v) + tt(u, u) 
= ir(u,u) — /iV r x a. 

□ 



Remark 4.1. Similarly, if u(t,x) is a smooth solution of the Euler equations (jl.4[) — (jl.5[) 
in Qj- satisfying the kinematic condition (|1.6j) . then the pressure p is determined by the 
Neumann problem of the Poisson equation: 

Ap = -V ■ (u ■ Vu) , d u p\ r = it(u,u). (4.4) 

4.2. Reduction to an IBVP with absolute boundary conditions via the Stokes 
equations. We will employ the solution of the unsteady Stokes equations to reduce 
our initial-boundary value problem into an initial-boundary value problem with abso- 
lute boundary conditions. More precisely, let u be a solution of the initial-boundary value 
problem of the Navier-Stokes equations (II. ip with nonhomogeneous boundary conditions 
(ll.6p - (|1.7p and initial condition (II. 3D . Let w be the solution of the nonhomogeneous 
initial-boundary problem for the unsteady Stokes equations: 

dtw = fiAw — Vq, V • w = 0, 

u;-L| r = 0, (V x w) 11 = a, (4.5) 
w\ t =o = u . 

Then v = u — w solves the following homogeneous initial-boundary value problem: 

dtv + (v + w) ■ V (v + w) = pAv — Vp, (4.6) 

V • v = 0, (4.7) 

v\t=o = 0, (4.8) 
subject to the absolute boundary condition: 

v ± = 0, (Vxv^O. (4.9) 

Therefore, if the initial-boundary value problem (14. 5p is well-posed with some appro- 
priate estimates of its solution, then the well-posedness for the nonhomogeneous initial- 
boundary value problem (jl.lj) - (jl.3j) and (ll.6p - (ll.7p is equivalent to the well-posedness of 
the homogeneous initial-boundary value problem (I4.6p - ()4.9p . 



4.3. Well-posedness of the Stokes equations and L 2 -estimates for the solutions. 

The initial-boundary value problem (|4.5p can be solved as follows. Taking the divergence 
in (|4.5I) . we can decouple the scaler function q by solving the Neumann boundary problem: 

Aq = 0, d u q\ r = -pV r x a, (4.10) 
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where the normal derivative of q follows from the kinematic condition in (|4.5p and V-w = 0, 
which implies from Lemma 12.21 that 

duQ\r = fJ,(Au,v) = — /A7 r x(Vxio)" = -/iV r x a. 

Since f T V r xadS = 0, there exists a unique solution q modulo a constant. We renormalize 
q so that J n qdx = 0. Therefore, it suffices to solve the linear parabolic equations: 

dtw — fj,Aw = —Vq, 

w x \ r = 0, (V x w) ]l =a, (4.11) 
w\ t =o = u , 

where q is given by (|4,10p . which is well-posed. In fact, by taking any smooth vector field 
A such that A ± \ T , = and (V x -A)"| r = a, and setting v := w — A, then the linear 
parabolic problem can be written as 

§- t v = ixAv + f , 
v ± = 0, (Vx«)"=0, 

where / := fiAA — Vq — ^mA, which is well-posed; its proof directly follows the arguments 
in Section 5. Furthermore, all the required estimates on w and q are available. 

We now make some necessary estimates for the solution uu of (|4.5|) . 
Taking the curl operation Vx in the equations in (|4.11|) . with g = V x w and h = 
V x g = -Aw, 

d t g = nAg; d t h = fj,Ah, 
together with the boundary conditions: 

5 ll |r = a, (V x g) 1 | r = V r x a, 

and 

h ± \ r = V r xa, (V x li) 1 | r = d t a. 

I 1 

We now make the L 2 -estimates for g and h. First note the following identities: 
d t {\g\ 2 ) = 2v(g, Ag), d t (\h\ 2 ) = 2fx(h, Ah). 
We integrate over Q to obtain 

-\\g\\ 2 2 = -2/i J^(g,V x(V x g)) dx 

= -2ii / (V x g, V x g) dx + 2fi / ((V xg) 11 x a, v) dS, 
Jn Jr 

which implies 

113111(^+2^/ ||Vxs|||(s)ds = [|^|||(0)+2/i / f ((Vxg)^xa, u) dSds for < t < T. 
Jo " " Jo Jr 

(4.12) 
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Similarly, we have 



dt u 112 



-2/i / (ft,, V x (V x h))dx 
Jn 

-2fi / (V x h, V x h) dx - 2 / ((d t a) x /i, ^) dS, 
Jn Jr 



which implies 



l(t) + 2n [ ||V x /i|||(a)ds = \\h\\l(0) -2 / [ ((d t a) x h, u) dSds for < i < T. 

(4-13) 

We now consider two different cases. 

Case 1: <9 t a = and y^a E L 2 (IY). In this case, we obtain from (|4.13p that 



\l(t) + 2n [ ||V x h\\l(s)ds 
Jo 

Since V • g = V • h = 0, we conclude 



2(0) < ||w(0,-)||h 2 = [|«o ll^- 



where C > is independent of 
Then we have 



l|V<?||£ 2(n) (i) +m||V^||£ 2(Qt) < C for t E [0,71, 



(Q T ) 



(4.14) 



2/i 



((V x 5) II x a, v) dSds 



Jr 



< 



< 



< 



2 M (e £ jf |V 5 | 2 dSds + C e ||a|| 2 2(rT) ^ 
M J J I V x 5 | 2 dx + C/x|| V 2 5|||2 ( n T ) + HHI|2(n T ) 

W / l v x 9\ 2 dx + C\\uo\\ 2 H 2 (n) + /i||a|| 
Jo Jn 



\h(r T )i 



where we have used the interpolation inequality in the second inequality above. 
Substitution this into (|4.12p yields 

WdWh^) + v\\Vg\\h { n t ) < c \\ u o\\ 2 H 2 {n) + Ml|a|||2 ( n T )- 
Combining (|4. 14|) ~ (|4. 15 j) with the fact V ■ w = V ■ g = V ■ h = 0, we have 

lkllfl-2 (n) (t) + M|h||| 3( n T ) <C , ||«oll|a(n) + ^lkllza(n T ) for 1 G t '^' 
where C > is independent of /i. 



(4.15) 



(4.16) 
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Case 2: dta 7^ and (a, dta) G L 2 {Tt)- In this case, we have 

/ {{dta) x h, v) dSds 
JT 



<e J j^\h\ 2 dSds + C £ \\d t a\\ 2 L2{TT) 

<nf [ |V x h\ 2 dxds + - [ [ \h\ 2 dxds + C\\d t a\\ 2 L2(n v 
Jo Jn V Jo Jn T 



where we have used V • h = and the interpolation inequality in the second inequality 
above. 

Substitution this into (|4.13p yields 

ft r< r-t 



l{s)ds. 



r c r 

\\h\\t(t) + H / || V x h\\ 2 2 {s)ds < \\u \\ 2 H 2 in ) + C\\d t a\\ L 2 {VT) + -± 

Jo A* Jo 

Then the Gronwall inequality yields 

\\hg{t)<M{v,T){\\u f H2m + \\d t a\\ 

i 2 (r T ))- 

Combining this with (|4.12|) . we conclude 

IHlif (n)(t) + W w W 2 HHn T ) < M{fx,T){\\u \\ 2 H 2 m + \\{a,d t a)\\ 2 L2{rT) ), (4.17) 
where M{fx, T) > depends only on \x and T. 

Proposition 4.1. Le£ w; G i/ 3 (f2) 6e a solution of problem (|4.5p , 27ien 

(i) If dta = and */7Ja G -L 2 (Tr) ; i/ien £/iere exists C > 0, independent of fi, such that 
(gJSH ZioZds; 

(ii) J/<%a 7^ and (a,d t a) G L 2 (T T ), i/tere exists M = M(p,T) > snc/i i/tai (|4TT7|) 
holds. 

5. A LINEAR INITIAL-BOUNDARY VALUE PROBLEM 

In this section, we deal with the absolute boundary problem of a linearized parabolic 
system. Let 0{t,x),O < t < T, be a smooth vector field on Q satisfying the absolute 
boundary condition: 



/i J 



r = (V x 0)11 

Define the pressure function p»{t,x) for < t < T by solving the Poisson equation: 



= 0. 

r 



f A Pf3 = -V-{{(3 + w)-X7(f3 + w)), 
\ 9 v pp\ r = iv{$ + w,/3 + w), 



subject to the normalization f^ppdx = 0. Here and hereafter, the time-variable t is 
omitted for simplicity in the equations if no confusion may arise. Thanks to Lemma 12. 11 



V • ((/? + w) ■ V{(5 + w)) dx = {((3 + w)- V(/3 + w),v)dS = - / n (/? + w, (3 + w) dS, 
n Vr Vr 

hence the above boundary value problem (|5,ip has a unique solution. The elliptic estimates 
yield that pp is smooth both in t and x. 
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We now consider the following initial-boundary value problem in Q: 

d t v + (j3 + w)- V(v + w) = Av - Vpp, (5.2) 
subject to the initial-boundary conditions: 

v\t=o = 0, (5.3) 



= LO' 

r 



= 0. (5.4) 

r 

In this section and Section 6, we denote uj = V x v again without confusion and always 
set [i = 1 without loss of generality for the existence proof. 

To solve problem f)5.2() — fj5.4[) . we recall that the Laplacian A acting on vector fields in 

D(A) = {v£ H 2 {n-R 3 ) : v ± \ r = J\ r = 0} 

is self-adjoint on L 2 (Q) and is negative-definite: 



(Av, v) = - / \V-v\ 2 dx- / \oj\ 2 dx < 
Jn Jn 



for every v £ D(A). The Laplacian with domain D(A) will be still denoted by A for 
simplicity. According to Hille-Yosida's theorem, A is the infinitesimal generator of a 
strongly continuous semigroup of contractions on L 2 (J7), denoted by ( e_tA ) t>0 - Indeed, 
A is the infinitesimal generator of an analytic semigroup, and the L 2 -domain -D(A) of A 
is invariant under e~* A . In addition, for each t, e _tA commutes with the curl operation 
Vx and the divergence operator V- (see [9] for the details). 

For every / £ D(A), tp(t) = e~ tA f solves the following evolution equation: 

d t ip = Aip, <p(0,-) = f. 
While, the regularity theory for parabolic equations yields that, if / is smooth, so is <p. 
Our next aim is to solve the homogenous equation: 

dt<p=(A-e-v)<p, v\ t =o = f, 

where / £ D(A), and 9 is a smooth vector field on £2 (independent of t). To this end, 
we show that A = A — 9 ■ V with domain D(A) = D(A) is the generator of a strongly- 
continuous semigroup. 

Lemma 5.1. Under the above notations, (A,D(A)) is a densely defined, closed operator, 
which is indeed the infinitesimal generator of an analytic semigroup on L 2 (Q). 

Proof. If a vector field (p satisfies the absolute boundary conditions, then 

\V<f\ 2 dx = - I A(\ip\ 2 )dx- I (A(p,(p)dx 

= \Jd v {\<p\ 2 )dS-j{Av,y)dx 



n((p,ip)dS— / (A(p,cp))dx 
r Jn 

< [ \ip\ 2 dS- [ (< A^,(p)dx 
Jr Jn 

< e\\Atp\\l + C\\<p\\l + \\A<p\\ 2 \\<p\\2, 
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that is, 

||v^||2<c(|M|2 + ||A^|| 2 |M| 2 ). 

Hence, we have 

11(0 • v^H 2 < ll^llLllv^ll 2 < cpwUMl + ||A^|| 2 |M| 2 ) < cwewlMt + £ ||a^|| 2 , 

for any e G (0, 1), where C(e, H^H^) is a constant depending only on e and H^H^. According 
to Kato's perturbation theorem (for example, see Theorem 2.1, page 80 in |38j), all the 
conclusions follow. □ 

The analytic semigroup with infinitesimal generator A — 9 ■ V with domain D(A) is 
denoted by (e* {AM - 9 ' v) ) t>0 . 

Next, we want to solve the nonhomogeneous evolution equation 

dtip= (A -{fi + w) -V)<p + h, 93(0,-)=/, 

where /3(t,x) is a smooth vector field and satisfies the absolute boundary conditions such 
that 

Mi)- P(s)\\oo < C\t - s\ a 
for some constant C and a G (0, 1], w(t,x) is the unique solution of (|4.5p . h(t, •) G L 2 (S7) 
for each t G [0,T], and / G D(A). 

Lemma 5.2. Let 



/)= sup ((A- (/? + «;)(*,-) ■V)^,ip) : iJGD(A M ) and ||^|| = 1 

0<t<T 

T/ien p < oo. If ip €z D(A) satisfies the Poisson equation: 

(A-(p + w)(t,-)-V-p-l)i/> = <p, t<T, 
and ip G L 2 (Q), then \\ip\\ < \\ip\\ and ||V^|| < C\\ip\\ for some constant C. 
Proof. It is easy to see that 

((A — (P + w)(t, •) • V)ip,4>) dx < + w)(t, OHoo / l^l 2 ^ 

so that p < oo. Since 

(ip,t/))dx = / ((- A + (0 + w)(t,-) • V + p + l)i>,ip)dx > 



which gives the first estimate in the lemma. Using the Bochner identity (|3.6p . we have 

\V^\ 2 dx = \ ! A(\iP\ 2 )dx- [ {Ai>,^)dx = \ I d v m 2 )dS- [ (AiP, iP)dx 
2 In .la 2 .lv .In 



ir(ip,ij>)dS- / (Atp,ip)dx 

\^\ 2 dS- I {<p+((p + w){t,-)-V + p+l)i/>,il>)dx 
Jn 

< e ||W||l + C|MH + ||^|| 2 ||V|| 2 + UP + w)(t, OlloollV^lb Mb + (p + 1)W1, 

which implies 

||v^|| <ch|. 



18 GUI-QIANG CHEN, DAN OSBORNE, AND ZHONGMIN QIAN 

□ 

Lemma 5.3. Let /3(t,x),0 <t<T, be a smooth vector field on Q satisfying the absolute 
boundary conditions (|1.14p . Suppose that 

+ IU <P for every t e [0,T] 

for some non-negative constant p. Let 

A(t) = A-((3 + w){t, •) • V - (p + 1)/ 

with domain D(A(t)) = D(A). Then, for each t, A{t) is the infinitesimal generator of the 
strongly continuous semigroup {e~( p+1 ' >s e s ^~(P +w ^ t '' r '^' ) } s >Q of contractions on L 2 (f2). 
Moreover, for every s,t, and r in [0,T], we have 

|| [A(t) - A(s)) ^(r)- 1 !! < C\\((3 + w)(t, •)-(/? + w)(s, -)||oo 
for some constant C . 

Proof. Let ip = A(r)~ 1 ip. Then tp G D(A) and, according to Lemma 15.21 ||V^|| < C||</?|| 
for some constant C. Since 

(A(t) - A(a))A(r)-V = ((/3 + w)(t, ■)-(/? + w)(s, •)) ■ W 

so that 

||(^)-A( S ))A(r)-Vl| < [|C0 + «j)(t,-)-(^ + «;)(s,-)||oo[|V^|| 

< C||(/3 + u;)(t,-) -(/? + «;)(«, Ollooll^l. 

□ 

Let us retain the notations as above. Assume that 

\\(f3 + w)(t,-)-(p + w)(s,-)\\ OD <C\t-s\ a for alls, t€ [0,T]. 
We want to solve the nonhomogeneous evolution equation: 

•) + (/? + «>)(*, •) • V^(t, •) = A^(t, •) + 

where ^(i, •) satisfies the absolute boundary conditions (|1.14j> . The above evolution equa- 
tion may be rewritten as 

d t V(t) = A(t)V(t) + e-( p+1)t h(t), 

where A(t) is given in Lemma 15.31 tp(t) = e^ p+1 ^^/(t), and ^(t) satisfies the absolute 
boundary condition. According to Theorem 6.1 in [38], page 150, there is a unique evolu- 
tion system U(t, s)onO<s<t<T such that the mild solution ^(t) is given by 

(i) = U(t, 0)^(0, •) + f e-( p+1)s U{t, s)h(s, -)ds 
Jo 

so that 

^(t, •) = e {p+1)t U{t, 0)^(0, •) + e {p+l)t f e-( p+1 >U(t, s)h(s, -)ds. 

Jo 

Therefore, we have the following theorem. 
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Theorem 5.1. Given the initial vector field f G D(A), h(t,-) G L 2 (£2), and f3(t,-) for 
t G [0, T] such that 

(3{t, •) G D(A M ) for each t G [0, T], 

supo<t<rll^(*>-)lloo < oo, 

+ w)(t,-)-((3 + w)(s,-)\\ 00 <C\t-s\ a V S ,iG[0,r], 



there exists a unique solution tjj which solves the absolute boundary problem on of the 
linear parabolic equation: 

dt*p + (/3 + w) ■ Vi> = Aip + h. 

6. Construction of local solutions 

The main goal of this section is to prove the existence of a unique strong solution, local 
in time, to the Navier-Stokes equations with the nonhomogeneous boundary conditions. 
To this end, we establish apriori estimate for solutions of the linear parabolic equations 

In order to state our results, we use the following norm: For a vector field v(t,x),0 < 

t < t,x g n, 



I«(V)IIn = ^Ht,-)\\ 2 H 2 + \\vt(t,-)\\ 2 m - 



6.1. Main estimate. Let T > be a fixed but arbitrary constant. Let f3(t, x), < t < T, 
be a given smooth vector field satisfying the absolute boundary conditions (|1.14p . Let 
v = V(J3) is the unique solution to the linear parabolic equations (|5.2|) - (|5.4p . with pp is 
the unique solution to problem (|5.ip and w(t, x) the unique solution of the initial-boundary 
value problem (|4.5p for the Stokes equations. Then the main apriori estimate is given in 
the following theorem. 

Theorem 6.1. Let 0(t,x),O < t < T, be a vector field which satisfies the following 
conditions: 

(i) For every 0<t<T, (3(t, •) G H 2 (n) and d t (3(t, •) G 

(ii) For each < t < T, (3 ± \ r = (V x 0)W\ T = 0; 

(iii) t —¥ \\P(t,-)\\ 2 H 2 + \\dt/3(t, -)\\ 2 H i is continuous; 

(iv) 0(0,0 =0. 

Let v = V(f3) defined by (|5.ip - (|5.4p . Then there exists a constant C depending only on 
the domain Q such that the following inequality holds: 

Ht,-)f N <C||u ||l 2 e^) + C r e W)-W)(l+ \m W )(s,.)\&) 2 ds, (6.1) 

J o 

where 

Q(t)= f {l + \mw)(s,-)\\ 2 N )ds. (6.2) 
J o 

To establish estimate (|6.ip , we will frequently use an elementary L 2 -estimate for elliptic 
equations, which can be stated as follows: If <f> is a solution to the Neumann boundary 
problem: 

= V • /, 



Mr = C/», (6 ' 3) 
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then we have the Solonnikov estimate: 

l|V#a < ||/||2. (6.4) 

Estimate (|6.4|) is a special case of Theorem 2.2 in Solonnikov [33]. In fact, estimate 
(|6.4|) follows easily from an integration by parts argument. Since 



IIV0H1 



f (cf), A(f>) dx + / 4>d v 4>dS = - [ ((f>,V-w)dx+ [ 0d u (pdS 
Jn Jr Jn Jr 

(V<f>,w)dx- [ <j)(w,v)dS+ [ <pd u (f>dS= [ (V</>,w)dS, 



so that (|6.4p follows from the Schwartz inequality directly. 

We now present the proof of Theorem l6.11 Throughout the proof, we will use C to denote 
a constant depending only on the domain Q, which may be different at each occurrence. 
We also omit the lowerscript f3 for simplicity: For example, we write p for pp. 

Let (1 = V't),w = Vxd, and ip = V x u for simplicity. Let ft denote the time-derivative 
dtf for a vector field /. Observe that vt again satisfies the absolute boundary conditions 
(|1.14p . Thus, according to Corollary 13.11 

\\v t \\ H i<C\\(v t ,d t ,u t )\\ 2 , \\v\\ H 2 <C\\(v,Vd,^)\\ 2 . (6.5) 

It is therefore natural to bound the time-derivative of each term appearing on the right- 
hand sides of (|6.5|) . However, as a matter of fact, we are unable to bound ^||Vfi|||. A 
different approach to handle ||Vd||2 is required. 

Taking the divergence V- on both sides of the linear parabolic equations (|5.2|) . one 
obtains 

d t d = Ad + V ■((f3 + w)-V((3-v)). (6.6) 

From (I5.2I)-(I5.4I). we have 



(Av, V) = — 7r(/3 + W, V + w) + 8 V p/3 = 7T (f3 + w, (3 - V ) 

so that, according to Lemma 12.21 d satisfies the Neumann boundary condition: 

d u d\ r = Tr(p + w,(3-v), (6.7) 

which is nonhomogeneous. We define, for every < t < T, a function q(t, •) by solving 
the Poisson equation: 

f Aq = -V-((P + w)-V(P-v)), 
\ 9 v q\ T = 7T (/3 + to, {} - v) , 

subject to j n qdx = 0. 

Let g = d — q. Then g solves the linear parabolic equation: 



(6.8) 



dt9 = Ag - q t , 
d u g\ T = 0. 

Since ||Vd|| 2 < \\(Vg,Vq)\\ 2 , we have 

||t;|| H2 < C\\(v,Vg,Vq,-, 
We thus consider the following function: 

F(t) = \\(v,^)(t r )\\l + \\(Vg,Vq)(t,-)\\l + \\(v t ,d t ,u t )(t, 



(6.9) 



(2- 
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Then |KV)Hn < F(t), and 

F(0) = ||(«o • V)u + VP0II2 + ||V x (n • V)«o||l, 
where po is the solution to 

Ap = -V • ((u • V)u ) , d v p \ r = ir(u , u ). 

Hence 

F(o)<c|Kl|^ 2 , 

and ()6,ip follows from the following lemma. 

Lemma 6.1. Under the assumptions in Theorem \6.1\ and notations introduced above, we 
have the following differential inequality: 

jF<C (||03, w)f N + 1) (F + 1 + ||03, ii;)^) • (6.10) 

Proof. Using (|5,2p - (|5,4p and integration by parts, we have 

= 2 f (v,Av)dx-2 f (v,Vp)dx-2 f (v, (f3 + w) ■ V(v + w)) dx 



dt 



n 



< 2||u||^ 2 +2||u|| 2 ||Vp||2 + 2||u|| 2 ||/? + w;|| 0O ||V(u + «;)||2. 
Using the Sobolev embedding, together with the L 2 -estimate for p, we have 

^H! < CII/3 + w\\ H 2 (F + 1) \\Vw\\ 2 . (6.11) 

Next we show the following inequality: 

^IKII 2 < C (1 + H/3 + HIn) (F + 1 + KUl + \\Vwg) • (6.12) 



By differentiating (15, 2p . we obtain 

d t v t = Av t -(/3 + w)- V(v t + w t ) - (fit + w t ) ■ V(v +w)- V Pt , (6.13) 

and vt also satisfies the absolute boundary condition (|1.14p . Since pt is the time-derivative 
of p, it solves the Neumann problem: 

A Pt = -V-{(p t + w t )-V(P + w))-V-{(P + w)-V(Pt + w t )), 
d uPt \ r = 2ir(/3 t + w t ,(3 + w). { °- LV 

According to the Solonnikov estimate (|6.4p . 

||Vpt|| < \\((3 t + w t )-V((3 + w) + V-(((3 + w)-V(f3 t + w t ))\\ 2 

< \\Pt + wt\k\WP + w)\\ 4 + + ^HoollVA + Vw t \\ 2 

< C\\(3 t + w t \\ m \\f3 + w\\ H 2, (6.15) 

where the second inequality follows from the Holder inequality, and the last one follows 
from the Sobolev imbedding: ||T|| 4 < C||T|| H i and \\T\loo < C||T|| H 2 in M 3 . Using ([633]) 
yields 

4lHli = 2 / {v t ,Av t )dx-2 [ (v t ,((3 + w)-V(v t + w t ))dx 
at Jn Jn 

-2 [ (v t ,V Pt ) dx-2 I (v t , (p t + w t ) • V(u + w)) dx. 
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Then, using the Holder inequality, we have 

J t W v t\\l < -2\\Vv t \\ 2 2 -2 J^ir(vt,v t )dS + \\(3 + w\\ 00 \\vt\\ 2 \\Vvt + Vw t \\ 2 

+2||u t || 2 || Vptlb + ||A + ^ilUll^tlUllVv + Vw\\ 2 , 

where we have used (|3.3p and (|3.4|) . According to the trace imbedding theorem (Theorem 
1.5.1.10 of Grisvard [12J, page 41), 



H 1 ' 



Thus we establish 



n(v t ,v t )dS\ < C\\v t u2 

Jr 

sh 

J t \\ v tf < C\\v t \\ 2 + 2\\v t \\ 2 \\Vp t \\ 2 + C\\p + w\\^ 2 \\vt + w t \\ 2 2 

+C\\p t + w t \\ H i \\v t \\ H i \\Vv + Vw\\ 2 , (6.16) 

which implies (|6.12p . 

Similarly, u> = V x v evolves according to the vorticity equation: 

d t iv = Alu-V x {{(3 + w) ■ V{v + w)), 
w"| r = 0, 

and ojt satisfies the evolution equation: 

f dtut = Auj t -V x {{p t + w t ) ■ V(v + w) + {f3 + w)- V{v t + w t )) , 

\ <4 

Moreover, V • u> t = 0. Hence, integration by parts yields 

d 2 
^lkl| 2 

= 2 (u t ,Au t )dx-2 / (wt,Vx {{(3 t + w t )-V{v + w) + {(3 + w)-V{v t + w t )))dx 
Jn Jn 

2 



= 0. 

r 



(6.17) 



(6.18) 



= -2 / | V x utf dx-2 / (V x uj u (fit + wt) ■ V(u + w) + (fi + w) ■ V(v t + w t )) dx 
Jn Jn 

< g II (A + «*) • V ( u + w) + G9 + w) • V(«t + u*)lll, (6.19) 



which yields 



^Ikiil < ci^ + hIn [f + IKv^v^nl) . (6.20) 



Next we consider HdtHl to show that 



±\\ dt \\ 2 < C\\P + w\\ 2 N {F + \\P\\ 2 N ) . (6.21) 



Indeed, differentiating equation (|6.6j) for d to obtain 

fltdt = + V ■ ((A + wt) ■ V{(3 -v) + (0 + w)- V{(3 t - v t )) , (6.22) 
subject to the Neumann boundary condition: 

d v <k\ T = n(fit + Wt,/3-v) + tt(P + *>,&- v t ). (6.23) 
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It follows that 

-||d t ||! = 2 / {d t ,Ad t )dx + 2 [ (d u V ■m + w t )-V((3-v) + (p + w)-V(Pt-v t )))dx. 
at Jn Jn 

Performing integrating by parts for the two integrals on the right-hand side and noting 
that the boundary integrals cancel out, we have 

^-\\d t \\ 2 = -2 [ \Vd t \ 2 dx + 2 [ (Vdt,(0 t + wt) • V((3-v))dx 
dt Jn Jn 



+2 / {Vd t ,(f3 + w)-V(Pt-v t ))dx 

< C\\(f3 t + w t ) ■ V(/3 - v)\\ 2 + ||09 + ti7)- V(A - v t )\\ 2 2 

< C||A + w t \\l\\V(3 - Vv\\l + + w\\l\\V(3 t - Vv t g 

< C||/3 + w\&[\\l3\\ 2 H 2 + \\v\\ 2 H 2 + ||VA - Vv t \\ 2 2 ), 

and (IOTP follows. 

tno corAnn.nrn or r\ ornra f nro Tn *a +" ic 1*70 n oorl r\r~ii inn 

dt 



Next we handle the second-order derivative. That is, we need to bound 4:\\Vg{t, -)|||, 



f t \\Vq(t,-)\\l, and f t U(t, -)\\ 2 L2 . We handle them one by one. 



Since g satisfies (|6.9p . integration by parts yields 

4||V<7||| = -2 / (Ag) 2 dx + 2 [ (Vq t ,Vg) dx < \\Vq t \\ 2 \\Vg\\ 2 . 
at Jn Jn 

While, qt solves the Poisson equation: 

Aq t = -V • ((A + wt) ■ V(/3 -v) + {(3 + w)- V(A - v t )) , 

subject to 

d v qt\ T = ir(f3 t + w t ,f3-v) + ir(f3 + w,(3t-vt) 
= -(V-((P + w)-V((3-v)),v). 

Thus, according to the Solonnikov estimate (|6.4p . 

||V?t||2 < \m + w t )-V(f3-u) + ((3 + w)-V((3t-v t )\\ 2 

< \\P + w\\k(\\P\\hi + Nla + l|V(A - vt)\\ 2 ) 

< \\p + w\\ N (\\/3\\ N + VF) , (6.24) 



and hence 



^[|V^(t,-)||l<2[|/3 + HlN(ll ?||n + v.F) v'.F. 



To estimate g|[|Vg(£, Olll; we begin with 

■||V?(t,-)||l=2 / (Vq t ,Vq)dx = -2 [ qAq t dx + 2 [ qd v q t dS. 

Jn Jn Jr 



d_ 

dt 1 
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Using the boundary condition for qt and integrating by parts again, we obtain 
d 
~dt 



i-||Vg(t,-)||| = 2 I qV ■((p t + w t )-V(p-v) + (p + w)-V(p t -v t ))dx 



+2 9 (tt(/% + w t ,(3-v)+ Tt(f3 + w,(3 t - vt)) dS 
-2 / (Vg, (fit + vt) • V(/3 - «) + 09 + w) ■ V(A - u t )) dS. 



Hence, we have 

d 



dt \\Vq\\l <C\\(3 + w\\ n {\\P\\ n + VF)Vf. (6.25) 

Finally, we establish a differential inequality for HVill- Recall that the vorticity u evolves 

according to the parabolic equation (|6.17p . Taking the curl operation Vx both sides of 
equation (|6.17j) . we obtain the evolution equation for ip: 

d t ip = Atp - V x V x ((/3 + w) ■ V(u + w)) , (6.26) 

so that 

d t (\iP\ 2 ) = 2(Ail>, ip) - 2(V x V x ((/3 + to) • V(u + w)) , 
Integrating over Q and performing integration by parts, one then obtains 
d 



jm,-)\\l = 2 / (A^)dx-2 / (VxVx {(p + w)-V(v + wj), i>) dx 
<ii Jn Jn 

V x ip\ 2 dx-2 / (V x (((3 + w) • V(v + w)), V x ip) dx 



+2 y x (V x V), i/) ^ - 2 y (V x (03 + w) ■ V(v + w)) x ^, dS, 

(6.27) 

where we have used the fact that V • ip = 0. Now we have to handle the last two boundary 
integrals. The vector identity: 

VxV = Vx(Vxtj) = -Auj, 

yields 

(ip x (V x V0> f) = (Aw Xip,v) = ((Aw)" X ip,u). 
However, since w"| r = 0, it follows from the vorticity equation (|6.17p that 

(Aw) 11 = (V x ((J3 + w)-V(v + w)))K (6.28) 
Therefore, the two boundary integrals sum up to zero. Hence, we have 

^-\\tp\\l = -2 [ |v x ip\ 2 dx-2 / (V x ({(3 + w) ■ V(v + w)),V x ip)dx. (6.29) 
Jn Jn 

To estimate the second integral, we consider a = V x (X • VY"), where X and K are two 
vector fields. Since 

(X ■ VY) k = X i V i Y k , 



then 
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a J = ie iW (v k (X a X7 a Y l ) - V l {X a V a Y k 

= e jkl (V k X a )(V a Y l ) + ^s ]kl X a (V k V a Y l - V,V a F fc ) 
= e jkl (V k X a )(V a Y l ) + X a V a (V x VJ . 



V x (X ■ VY) = e jkl (V k X a )(V a Y l ) + X ■ V (V x Y) 



That is, 
so that 

|V x (X • VY)\ < |VX||VF| + |X||V (V x Y) \. 
Using this inequality in (|6.29p . one deduces that 

— < -2 f |V x ^\ 2 dx + 2 [ |V x VI \V(J3 + w)\ |V(v + w)|dx 



+2 / | V x VI |/? + to | \V(oj + V xw)\dx 
Jn 

< -2||V x V||i + 2||V x V||2||V(/3 + u»)|| 4 ||V(v + u»)||4 
+2||V x Vl|2||/3 + ^||oo||Vw + V(V x w)\\ 2 

< -2||V x ij\\l + C\\V x t/j\\ 2 \\P + w\\ H 2\\v + w\\ H 2, 

'' "\l*<C\\l3 + w\\^{F+\\w\\l 2 ). (6.30) 



which implies that 

dt 1 

Combining all the estimates yields (|6.10p . We thus complete the proof of Theorem 
EH □ 



6.2. Local strong solutions. For T > 0, denote Wt the space of all vector fields /3(t, x) 
that satisfy conditions (i)-(iii) in Theorem 16.11 with the uniform norm: 



w T = 

0<t<T 



sup v / II/3(*.0II^ + IIW,0IIhi- 



Then is a Banach space. According to (|6.ip . there is a constant C > depending 
only on Q such that 



\\V(P)\\w T < C^Coe 

+ CVTe C{1+mw) ^ T (l + 11(^^)11^). (6 . 31) 

Equipped with the apriori estimate (|6.ip . we are now in a position to establish the ex- 
istence of a local (in time) strong solutions for the initial-boundary value problem with 
nonhomogeneous boundary conditions. 

Theorem 6.2. Given Cq > 1, there exists T > depending only on Co, the viscosity 
coefficient \i > 0, and the domain £2 such that, for any given Uq £ H 2 (£l) which satisfies the 
absolute boundary conditions (ll.6p - (|1.7[) . ||uo||ff 2 < Co> andV-uo = 0, there exists a strong 
solution u(t, x) of (]l.ip - (ll.3p and (|1.6|l - (jl.7|) with the form u(t, x) = v(t, x) +w(t, x), < 
t <T, such that 

(i) For every 0<t<T, v(t, •) G H 2 (tt) and d t v(t, •) 6 H 1 ^); 
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(ii) For every < t <T, v L \ r = and (V x u)H | = 0; 

(iii) t — ► \\v(t, -)\\ff2 + ||5fc u (^?')[|jfi *s continuous; 

(iv) f |j=o = 0, and w satisfies the Navier-Stokes equations 

dtv + (v + w ) • V(v + w) = Av — Vp, V • v = 0, 

where p solves the Poisson equation for each t £ (0, T]. - 

Ap = —V • {{v + w) ■ V{v + w)) , d u p\ T = ir(v + w,v + w), 

and w solves the initial-boundary value problem (|4.5f) for the unsteady Stokes equations. 

Proof. Step 1. If 0i(t,x) and (3 2 (t,x), < t < T, are two vector fields which satisfy 
conditions (i)-(iv) in Theorem 16.11 then U := v\ — V2 with = V{j3k),k = 1,2, satisfies 
the linear parabolic equations: 

d t u = AU - (A + w) ■ VU — VP — (A - 2 ) • V(v 2 + w), U\ t=0 = 0, 
subject to U ± \ r = and (V x {/)" = 0, where P = p^ — pp 2 is the unique solution to 
AP = -V • ((/3i + w) • V(/3i + w) - (p 2 + w) • V(/3 2 + it;)) 

such that 

C^Plp = 7T (/3l + W, fix + W) — 7T ((3 2 + W, f3 2 + w) . 

Applying the same arguments used in the proof of Theorem 16. 1\ we obtain that there 
exists C > depending only on the domain $7 such that the following inequality holds: 

\\V(Pi)-V(/3 2 )\\w T 

<C^e C{l+WlM ^^ (6.32) 



It follows from (|6.3ip - (|6.32p that, for every K > ^JCCq, there exists T > depending 
only K, C, and ||w||n such that, if v G Wj- with \\v \\w T — then ||V(v)||w T — an d 

2 

||V(ui) - V(«i)||w T ^ 3F1 - u 2||w T 

for any v\,v 2 G Wf such that t>i(0, x) = i>2(0, x) = 0, ||uj||w T < K , i = 1, 2. Therefore, 
there is a unique fixed point v £ Wy such that V(v ) = v. Then u = v + w, according to 
Theorem l6.H a solution of the nonhomogeneous initial-boundary value problem (|l.ip - (|1.3j) 
and (|1.6p - (|1.7p for the Navier-Stokes equations. 

Step 2. Incompressibility: If v(t, -),0 < t < T, is a fixed point of the velocity map V, 
then v is a solution to 

V ■« = (). (6.33) 
Setting /3 = i> in (|5.ip - (|5.2p and taking divergence to obtain 

<%(V • v) + V • ((v + w) ■ V(v + w)) = A(V • w) - Ap,, 

so that d = V ■ v solves the heat equation: 

8 t d = Ad, d(0, •) = 0. 

According to our assumptions: 

d t v + (v + w) ■ V (v + w) = Av - Vp v . (6.34) 
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Identifying the normal part of each term of equation (|6.34p . together with the boundary 
condition: v ± \ r = w L = 0, we conclude 

(An) 1 = {Vpv^ + ((v + w)-V(v + w)) 1 = 0. 

On the other hand, by Lemma 12.21 together with the boundary condition: w"| r = 0, we 
have 

d v d\ T = (Av,v) = 0. 

By the uniqueness of the Neumann problem, d(t, •) = for all t. □ 

7. Inviscid Limit in $7 c R n ,n > 3 

In this section we study the inviscid limit of the solutions to the nonhomogeneous initial- 
boundary value problem (|l.ip - (|1.3j) and (|1.6p - (|1.7p for the Navier-Stokes equations. 

Let n^ be the solution to the initial-boundary value problems (|l.ip - (jl.3p and (|1.6|) - 
(jl.7p for the Navier-Stokes equations for each fj, > 0. Let n be the solution to the initial- 
boundary value problem (jl.4p - ([1.6p for the Euler equations. Notice that all solutions n M 
subject to the same boundary conditions: uj[\ r = and (V x ii M )"| r = V r x a for the 
given smooth vector field a, while the solution n satisfies only the kinematic condition 
u ± \ r = which is independent of the viscosity constant u. 

Theorem 7.1. Let a £ L 2 ([0,T]; L 2 (T)) be a smooth vector field on T. Suppose that, for 
all [i G (0, no], a unique strong solution of problem (|l.ip - (jl.3p and (|1.6p - (|1.7p and a 
unique strong solution u € H 2 {^1) of problem (|1.4p ~( fl~6j) exist up to time T* > 0. Then 
there exists C > depending on fj,Q, T, |H|z,2(r ,T];£ 2 (r)j and \\u\\H 2 r\W 1 >°°(n)> independent 
of /J,, such that, for any T G [0, T*], 

sup \\u„(t,-)-u(t,-)\\ 2 <C(T)^^0 asnjO, (7.1) 

0<t<T 

and 

\\V(u u -u) (s,-)\\lds < C. 





It follows that the solutions of problem (|l.ip ~ (|1.3p and (|1.6p ~ (jl.7p for the Navier-Stokes 
equations converge to the unique solution u(t,x) of problem (fLUl - flTUp in L°°([0, T];L 2 (Q)). 

Proof. Let v n = — u. Then Vn satisfies the following equations: 

dtv^ = fj-Avf, - (Vf, + u) ■ Vv^ - VP M -v fl -Vu + fiAu, , . 

v-^ = o, [L2) 

and the initial condition: 

u„(0,.) = 0, (7.3) 

where = — p. 

Since both n M and u satisfy the kinematic condition (jl.6p . so does v^. Thus, by means 
of the energy method, we obtain 



d 



"77 H ' /< ih 2/i / v^-Av^dx- I {v^ + u)-V{\v^\ 2 )dx-2 I X P fl . < 
at Jn Jn Jn 

-2 / {v^ ■ Vn, Vfj) dx + 2u / (An, v M ) dx. 

Jn 
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Integrating by parts in the first three integrals, one then deduces that 
d 



= -2/x / |V x v fl \ 2 dx + 2[i / (vp x (V x ,v)dS 
Jn Jr 

-2 / (v^ ■ Vu, Vfj) dx + 2/j, / (Au,v^)dx 
Jn Jn 

= -2n\\V Vlx \\l - 2n J n(y li ,v IM )dS + 2n J (v^xb,is)dS 

-2 / (v fJ/ -Vu,v IJ ,)dx + 2n / {Au,v^)dx, 
Jn Jn 

where b = a — V x u. Furthermore, we use the following estimate: 

j^v^ x b,v)dS < 11611^^)11^1^2^) < C(\\v^\\ 2 L 2 iT) + \\a- V x u||| 2(r) ) 
to obtain 

^IKIll < -2/x||V^||l + 2/xC (||u M |||2 (r) + \\a- V x Hli^r)) 

+2||V«||oo||umII2 + 2mI|A«||2||u m ||2 • (7.4) 
Finally, we use the Sobolev trace theorem: 

2C||v / j£ 2 (T) < ||Vi; M ||| + C , ||u j[1 ||2, ||V x u||| 2 (r) < IMI^^) 
to establish the differential inequality: 

^IKIll + /*l|Vv At ||l < C ((||V«||oo + /*o)Wll + K\\ a \\h(T) + IMI^W) ■ ( 7 - 5 ) 
The Gronwall inequality yields that, for t G [0, T*], 

\K(t,-)\\t < MC^ e C ^(ll V ^ll- + ^- s »(||6||| 2(r) + \Hs,.)f HHn) )ds := Cn. (7.6) 
Hence 

rt 

Vv^s,-)\\ 2 ds<C, 



I 

Jo 



10 

which imply the conclusions of the theorem. □ 

In order to ensure the convergence of to u in the strong sense (say, in H 2 (Q)) up 
to the boundary, a necessary condition is that u must match with the boundary data 
(V x u) ± \ r = a. 
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